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Abstract 


Floods are a common natural disaster that can endanger urban and rural infrastructure. The purpose of this 
study was to prepare and analyze the flood risk map of Khancheh and Baronagh villages along a part of 
Suk Cham river. First, the topographic map was created using Google Earth and Global mapper. Then, 
using a topographic map made in AutoCAD Civil 3D software, the surface was created for the study 
interval. In the next step, using the constructed surface of cross sections, longitudinal profiles, river 
channel, right, and left flood banks of Suk Cham river were constructed in AutoCAD Civil 3D software 
environment. To simulate floods at return periods of 25, 50, 100, and 200 years, the cross sections were 
transferred to HEC RAS software. The Phil-Steiner method was used to convert the maximum daily flow 
rate to the maximum instantaneous flow rate according to the P-value of the Kolmogorov-Smirnov test 
(0.9). The results of flood frequency analysis showed that Pearson type 3 log distribution is the most 
suitable distribution. Based on the results, the flooding level for the return periods of 25, 50, 100, 200 
years was 0.77, 0.9, 1.15, 1.47 km’, respectively. According to the flood risk map in the 100-year return 
period (high risk area), about 41.596 and 4.596 of the two villages of Khancheh and Baronagh had at 
moderate to very high risk of flood, respectively. Establishing coastal walls and organizing the river along 
the river route through these two villages is essential, which can mitigate the risks of floods. The use of 
prepared plans for flood control measures, flood forecasting and warning systems, flood insurance and in 
general for structural and non-structural flood control measures in the region is recommended. 


Keywords: Flood risk map, Khancheh and Baronagh villages, Hydraulic model, AutoCAD Civil 3D, Suk 
Cham River Watershed 


* , Corresponding author: Reza Ghazavi E-mail: ghazavi@kashanu.ac.ir Tel: + 989171887058 


How to cite this Article: Sayyad, D., Ghazavi, R., Omidvar, E. (2021). Preparation and analysis of flood 
risk map using HEC RAS and RAS MAPPER hydraulic model (Case study: Sok Cham river of Kashan). 
Journal of Geography and Environmental Hazards, 10(3), 19-37. 

001: 1.69554. 1038 

Journal of Geography and Environmental Hazards are fully compliant 

with open access mandates, by publishing its articles under Creative 

Commons Attribution 4.0 International License (CC BY 4.0). 


Creative Commons Attribution 4.0 International License (CC BY 4.0) 


Geography and Environmental Hazards 
Volume 10, Issue 3 - Number 39, Fall 2021 
https://geoeh.um.ac.ir 


5 https://dx.doi.org/10.22067/geoeh.2021.69554.1038 
۱۹-۳۷ جغرافیا و مخاطرات محیطی. سال دهم شمارة سی و نهم پاییز ۱۶۰ صص‎ 
مقاله پژوهشی‎ 
RAS MAPPER و‎ HEC RAS و تحلیل نقشه خطر سیلاب با استفاده از مدل هیدرولیکی‎ ag 
(مطالعه موردی: رودخانه سوک چم کاشان)‎ 
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ce 


سیل یکی از بلاهای طبیعی رایج است که می‌تواند زیرساخت‌های شهری و روستایی و جان و اموال را با 
خطر مواجه کند. هدف از انجام این پژوهش تهیه و تحلیل نقشه خطر سیل دو روستای خنچه و بارونق در 
امتداد بخشی از رودخانه سوک چم است. ابتدا نقشه توپوگرافی با استفاده از Google Earth‏ و 
Global mapper‏ ساخته شد. سپس با استفاده از نقشه توپوگرافی ساخته شده در محیط نرم‌افزار 
AutoCAD Civil 3D‏ برای بازه مطالعاتی سطح ایجاد گردید. سپس با استفاده از سطح ساخته شده 
مقاطع عرضی. پروفیل طولی. کانال رودخانه. سواحل سیل‌گیر راست و چپ رودخانه سوک چم در محیط 
نرم‌افزار AutoCAD Civil 3D‏ ساخته شد. برای شبیه‌سازی سیلاب در دوره بازگشت‌های ۲۵ Or‏ 
Yer ۰‏ ساله مقاطع عرضی ساخته شده به نرم‌افزار HEC RAS‏ انتقال داده شد. برای تبدیل دبی 
حداکثر روزانه به دبی حداکثر لحظه‌ای با توجه به مقدار P-value‏ آزمون کولموگروف-اسمیرنوف C/A)‏ 


از روش فیل - استینر استفاده شد. نتایج حاصل از تحلیل فراوانی سیلاب نشان داد توزیع لوگ پیرسون 


Email: ghazavi € kashanu.ac.ir ۰۹۱۷۱۸۸۷۰۵۸ نویسنده مسئول‎ ۱ 


نحوه ارجاع به این مقاله : 
صیاد, JU ls‏ قضاوی, رضا, امیدوار, ابراهیم. Ort)‏ تهیه و تحلیل نقشه خطر سیلاب با استفاده از مدل 
هی درولیکی HEC RAS‏ و RAS MAPPER‏ (مطالعه موردی: رودخان» سوک چم کاشان). جغرافیبا و 
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نوع Y‏ مناسب‌ترین توزیع است. بر اساس نتایج حاصل» سطح سیل‌گیری برای دوره باز گشت‌های ۵ Or‏ 
در دوره بازگشت ۱۰۰ ag) JU‏ پرخطر) به ترتیب حدود 1۱/۵ و 1/۵ درصد از دو روستای خنچه و 
بارونق در معرض خطر سیلاب متوسط تا خیلی زیاد قرار دارند. ایجاد دیواره‌های ساحلی و سمان‌دهی 
رودخانه در طول مسیر عبور رودخانه از این دو روستا امری ضروری است که می‌تواند ol h>‏ ناشی از 
سیل را تعدیل نماید. استفاده از نقشه‌های تهیه شده برای انجام اقدامات کنتصرل سیل سامانه‌های پیش‌بینی 
و هشدار سیل. بیمه سیل و به‌طورکلی جهت اقدامات سازه‌ای و غیرسازه‌ای مهار سيل در منطقه پيشنهاد 
می‌شود. 

کلیدواژه‌ها: نقشه خطر سیل. روستاهای خنجه وبارونق. مدل هیدرولیکی. AutoCAD Civil 3D‏ 


حوضه آبخیز رودخانه سوک چم. 
۱- مقدمه 


سیل به جریان‌های فراتر از ظرفیت رودخانه‌های طبیعی و يا کانالهای مصنوعی گفته می‌شود که به اراضی 
اطرافب گسترش ly‏ کرده ‏ می‌تواند ایجاد عسارث کند D‏ ۲۰۲۰):میل یکی از بلاهلی:طبیغی alb‏ است که 
می‌تواند تأثیر زیادی روی زیرساخت‌های شهری و روستایی و ساختمان‌ها و جان و اموال مردم داشته باشد و به‌طور 
غیرمستقيم روی اقتصاد کشور اثر بگذارد (ایبراهیم " و همکاران: ۲۰۱۶؛ کیم " و همکاران ۲۰۲۰). در سال‌های اخیر سیلاب 
۷ درصد از بلایای طبیعی را به خود اختصاص داده است که به یک فاجعه طبیعی در سطح جهان تبدیل‌شده است و 
Gol cl: le E ally‏ ۱۱/۱ مرس از لفات انسانی در تسطح od E oes‏ 
ورد 3/3 15 (T4513)‏ 

با توجه به پراکنش نامنظم بارندگی در اغلب نقاط کشون هرساله در نقاط مختلف کشور سیلاب‌های شدید و 
تا ای اد موی کف بات ae Seed‏ وان oo agb‏ کر ر ریق شیر deles‏ نی و 
پهنه‌بندی lam‏ سیل برای رودخانه‌ها؛ یکی از روش‌هایی است که می‌تواند به مدیریت سیلاب و کاهش ریسک 
مناطق سیل گیر کمک کند. مدل‌سازی سیل یکی از روش‌های ارزیابی خطرات سیل برای نمایش بزرگی سیل در دوره 
باز گشت‌های مختلف است fl‏ و همکاران (A‏ 
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سال دهم تهیه و تحلیل نقشه خطر سیلاب با استفاده از مدل .... Y!‏ 


لازمه تحلیل سیلاب در دوره بازگشت‌های مختلف اندازه‌گیری و يا برآورد میزان دبی‌های اوج لحظه‌ای است. 
ازآنجایی که در اکثر رودخانه‌های leper Ol ud‏ در حوضه‌های خشک و نیمهعشک برآورد دبی‌های حداکثر لحظه‌ای 
ممکن نیست و با توجه به نقش این دبی‌ها در پهنه‌بندی سیل در دوره بازگشت‌های مختلف» آمار دبی‌های حداکثر 
لحظه‌ای با استفاده از روش‌های غیرمستقیم (فرمول‌های تجربی) از میانگین روزانه دبی در این حوضه‌ها محاسبه می- 
شود (شعبانلو و همکاران,۱۳۹۱). استفاده از مدل‌های هیدرولیکی بستگی به هدف مدل‌سازی می‌تواند یک‌بعدی یا 
دوبعدی باشد ولی استفاده از مدل‌های هیدرولیکی یک‌بعدی معمولاً رایج‌تر است. در مدل‌سازی یک‌بعدی مقاطع 
عرضی در طول شبکه جریان ترسیم می‌گردد و شبیه‌سازی جریان در طول مسیر رودخانه انجام می‌شود (وجتک و 
همکاران (QM‏ 

مدل HEC RAS"‏ یکی از پرکاربردترین مدل‌های تجزیه‌وتحلیل = Ob‏ در رودخانه‌ها و کانال‌ها و ترسیم نقشه- 
های پهنه‌بندی سیل در دوره بازگشت‌های مختلف است (ختاک " و همکاران. Y YA‏ تا به امروز پژوهش‌های 
متعددی در ایران و جهان درزمینه پهنه‌بندی و حطر سیل صورت گرفته است. 

ازجمله این پژوهش‌ها می‌توان به پژوهش‌های زراعتکار و حسن‌پور (۱۳۹۵) در پژوهشی با عنوان شبیه‌سازی 
سیلاب شهری بیرجند با استفاده از مدل HEC RAS‏ و GIS‏ به این نتیجه دست یافته‌اند که ميزان دبی اوج سیلاب- 
های با دوره بازگشت ۲۵ و ۱۰۰ سال به ترتیب ۲۰۷/۱۳ 17۳/۵7 مترمکعب بر Ad‏ به‌دست‌آمده است. 

حجازی و همکاران (۱۳۹۸) در پژوهشی با عنوان پهنه‌بندی خطر وقوع سیلاب در حوزه آبریز ورکش چای با 
استفاده از مدل HEC RAS‏ و HEC-GeoRAS‏ به این نتایج دست‌یافتهاند که سیلاب‌های با دوره بازگشت ۲۵ و ۵۰ 
ساله به ترتیب ۳" و ۱۰۰ کیلومتر از مساحت کل حوزه آبریز را تحت پوشش قرار می‌دهند. 

قنواتی و همکاران (۱۳۹۸) در پژوهشی با عنوان شبیه‌سازی پهنه‌های سیل گیر در حوضه آبخیز کن در شهر 
تهران به این نتایج دست یافته‌اند که در بالادست پهنه‌های سیل گیر گسترش عرضی کم‌تر و دبی زبادتری (ارتفاع 
سطح آب) دارند اما در بخش‌های کم شیب دبی کم‌تر (ارتفاع سطح آب) ولی پهنه‌های سیل‌گیر گسترش عرضی 
dioses‏ 

مهرورز و همکاران (۱۳۹۹) در پژوهشی با عنوان شبیه‌سازی سیلاب‌های رودخانه دره ائورت با استفاده از مدل 
هیدرولیکی HEC RAS‏ در GIS lass‏ رودخانه ارس به این نتیجه دست یافته‌اند که حداکثر عرض پهنای سیل گیر 
برای سیلاب با دوره باز گشت ۲۵ سال حدود za ۸٩7‏ است و حدود ۱۷٤۳‏ هکتار از اراضی دشت سیلابی زیرآب 


] Vojtek 
2 Hydrologic Engineering Center's River Analysis System 
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فرو می‌رود. اشاره کرد. دورجی ‏ و همکاران (۲۰۱۹) در پژوهشی باهدف مدل‌سازی سیل رودخانه آموچو با استفاده از 
HEC RAS‏ و GIS‏ به این نتایج دست یافته‌اند که شهر فونتشو لینگ در برابر سیلاب‌ه ای با دوره بازگشت YO‏ ۱۰۰ 
ساله آسیب‌پذيرند. 

اونگدز" و همکاران (۲۰۲۰) در پژوهشی باهدف تولید نقشه‌های خطر سیل در قزاقستان بیان کردند که روستای 
ولگو در رویداد سیل ۱۰۰ ساله سیل‌زده شد. آینالم " و همکاران (۲۰۲۰) در پژوهشی در رودخانه موگا مساحت‌های 
سیل گیر برای دوره بازگشت‌های ۵ ۸۰ ۲۵ Oe‏ ۱۰۰ به ترتیب ۱۸ ۰۲۱0۲۱ ۳۶ EY‏ کیلومترمربع گزارش کردند. 

تهیه و تحلیل پهنه‌های حطر سیل با توجه به نقش آن در مدیریت سیلاب‌ها امری ضروری در جهت مدیریت 
سیلاب دشت‌ها است. همچنین نقشه‌های خحطر سیل برای اقدامات سازه‌ای و غیر سازه‌ای و کنترل سیل ضرورت 
دارد. هدف از انجام این پژوش تعبین پهنه‌های سیل OS‏ تولید نقشه خطر سیل و تحلیل خطر عناصر در معرض 
خطر سیل دو روستای خنچه و بارونق با استفاده از نرم‌افزارهای Global Auto CAD Civil 3D. HEC RAS‏ 


و RAS Mapper‏ است. 


-Y‏ مواد و روش 

۱-۲- منطقه موردمطالعه 

حوضه آبخیز رودخانه سوک چم با مساحت حدود ۲۱۵/۸۳ کیلومترمربع یکی از حوضه‌های مشرف به شهر 
کاشان است که در استان اصفهان و مرکزی ایران واقع شده است. این منطقه بین طول جغرافیایی ۵۱ درجه و ٩‏ دقیقه 
و ۵ ثانیه شرقی و عرض جغرافیایی ۳۳ درجه و ۵۸ دقیقه و ۳۷ انيه شمالی واقع شده است. رودخانه سوک چم 
مهم‌ترین رودخانه این حوضه بوده و مسر حرکت آن از بالادست به پایین‌دست شهر کاشان و سرانجام به دریاچه 
نمک در کویر مرکزی ایران ختم می‌شود. روستاهای خنچه و بارونق در مجاورت بخشی از این رودخانه els‏ شده 
است. همچنین طول بازه مطالعاتی در این پژوهش از بالادست به پایین‌دست حدود ۱۱ کیلومتر است. در شکل ۱ 
موقعیت منطقه و بازه مطالعاتی آمده است. 


1 Dorji 

2 Ongdas 

3 Aynalem 

4 civil engineering design software 
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شکل ۱- موقعیت منطقه و بازه مطالعاتی 
-Y-Y‏ روش تحقیق 
M Y-Y‏ - دبی‌های حداکثر لحظه‌ای 
ایستگاه هیدرومتری بازه موردمطالعه فقط دبی‌های حداکثر روزانه را قرائت می‌کند ولی برای شبیه‌سازی سیل با 
دوره بازگشت‌های مختلف در مدل HEC RAS‏ به دبی‌های حداکثر لحظه‌ای نیاز است. در این پژوهش برای برآورد 


دبی‌های حداکثر لحظه‌ای با استفاده از آمار دبی‌های حداکثر روزانه ایستگاه هیدرومتری خنچه از دو روش تجربی 
فیل - استینر و سنگال به ترتیب رابطه ۱ و É‏ استفاده شد (فیل و استینر ؛ ۲۰۰۳؛ سنگال " ۱۹۸۳). 


Qmax=(0.8Q,+0.25)Q,+Q;)) /k 6 (رابطه‎ 


1 Fill & Steiner 
2 Sangal 
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× فاکتور تصحیح است که از رابطه Y‏ و ۳ به دست می‌آید: 


K=0.9123 6 )۲ (رابطه‎ 
X=(Q,+Q,)/2Q, ( (رابطه‎ 
Qmax-(4Q,-Q, -Q.)/2 04 (رابطه‎ 


مترمکعب بر ثانیه) روزی که دبی حداکثر در آن رخ‌داده است و ,۵ و Q,‏ به ترتیب میانگین دبی روزانه (برحسب 


-Y-Y-Y‏ تحلیل فراوانی سیلاب 

آنالیز فراوانی سیلاب با توجه به مقادیر دبی‌های حداکثر لحظه‌ای برآورد شده از فرمول‌های تجربی فیل - استینر و 
سنگال با استفاده از ۵ توزیع آماری نرمال» لوگ نرمال» لوگ نرمال ۳ پارامتره مقادیر حد و لوگ پیرسون تیپ ۳ در 
نرم‌افزار EasyFit‏ برازش داده شد و با توجه به آزمون کلموگروف اسمیرنوف و مقادیر P-value‏ (در سطح اعتماد ۵ 
Go»‏ بهترین توزیع و روش تجربی انتخاب شد. سپس دبی با دوره بازگشت‌های ۵ Oe‏ ۰۱۰۰ ۰۲۰۰ ساله 
برحسب مترمکعب بر ثانیه با استفاده از بهترین توزیع آماری و داده دبی حداکثر لحظه‌ای حاصل روش تجربی برآورد 
NIE‏ 

۳-۲-۲- ساخت داده‌های هندسی 

ساخت داده‌های هندسی (ترسیم کانال رودخانه» پروفیل طولی رودخانه, سواحل دشت سیلابی چپ و راست و 
ترسیم مقطع عرضی و جهت جریان) اولین گام شبیه‌سازی سیلاب در دوره بازگشت‌های مختلف با استفاده از مدل 
HEC RAS‏ است. در این پژوهش برای ساخت داده‌های هندسی از نرم‌افزارهای» Global Mapper «Google Earth‏ 
AUTO CAD civil 3D‏ (نرم‌افزار مهندسی در ارتباط با طراحی مسیر) استفاده شد. در ابتدا با استفاده از نرمافزار 
google earth‏ و global mapper‏ نقشه توپو گرافی برای منطقه موردمطالعه ساخته شد. در شکل ۲ نقشه توپوگرافی و 
سطح ساخته شده در نرم‌افزار 30 AUTO CAD civil‏ نشان osla‏ شده است. سپس با انتقال نقشه توپوگرافی ساخته 
شده به نرمافزار AUTO CAD civil 3D‏ سطحی از منطقه موردمطالعه ساخته شد و از روی سطح ساخته شده 
رودخانه و کانال آن و مشخصات هندسی رودخانه (مسیر رودخانه جهت جریان. سواحل سیل گیر چپ و راست؛ 
مقاطع عرضی) ترسیم شد. سپس از داده‌های هندسی ترسیم شده در نرم‌افزار AUTO CAD civil 3D‏ خروجی 
گرفته و به نرم‌افزار HEC RAS‏ انتقال داده شد. در شکل ۳ داده‌های هندسی ترسیم شده (مسیر رودخانه. جهت 
جریان» سواحل سیل گیر چپ و راست. مقاطع عرضی) در نرم‌افزار AUTO CAD civil 3D‏ نشان داده شده است. 


سال دهم ang‏ و تحلیل نقشه daz‏ سیلاب با استفاده از مدل .... Yo‏ 


مزیت این نوع ترسیم از داده‌های هندسی با استفاده از نرمافزار AUTO CAD civil 3D‏ دقت GUL‏ این نرمافزار در 
ترسیم داده‌های هندسی Aud‏ 


3430314 
Mele 


Sr 18 Oe 


شکل Y‏ - ساخت نقشه توپوگرافی ساخته شده (سطح) با استفاده از Auto CAD Civil 3D‏ 


YA‏ جغرافیا و مخاطرات محیطی شمارة سی و نهم 


شکل ۳- الف) موقعیت دو روستا خنجه و بارونق ب و پ) داده‌های هند سی ایحاد شده در نرم‌افزار 
(c Auto CAD Civil 3D‏ ترسیم سه‌بعدی کانال 3 مقاطع عرضی در HEC RAS‏ 


سال دهم تهیه و تحلیل نقشه daz‏ سیلاب با استفاده از مدل .... yv‏ 


4-۲-۲- تعیین ضریب زبری مانینگ 
برای تعیین ضریب زبری مانینگ روش‌های مختلفی ازجمله قضاوت کارشناسی, بازدیدهای میدانی. جداول 
پیشنهادی بر اساس پژوهش‌های سایر محققین وجود دارد. در این پژوهش ضریب زبری مانینگ بر اساس جدول 


پیشنهادی چو جدول ۱ و بازدید صحرایی به دست آمد (چو À‏ 404( 


جدول ۱- جدول چو برای ضریب زبری مانینگ 


. ۳ ضریب زبری 
ردیف گروه توضیحات 
حداقل متوسط حداکثر 
۱ 0 ۳ ۰۳۳ تمیز و مستقیم و فاقد مخزن ذخیره 
J^ ۰/۳۵ JY Y‏ همانند شماره ۱با سنگ و علف هرز بیشتر 
کانال اصلی 7 
LEA ۰/۳۳ ۳‏ ۰۱/۰۹۵ تمیز و مارپیچ و بعضا دارای مخزن ذخیره 
۶ ۰/۳۵ ۰:۵ ۰/0 همانند شمار Y‏ همراه با علف هرز و سنگ 
0 ۵ ۳ 0 مرتع با علف‌های کوتاه 
lL‏ دشت 
“LYO ۳ 1‏ ۰/۵ مرتع با علف‌های بلند 


۵-۲-۲- اجرای مدل یک‌بعدی HEC RAS‏ در شرایط ماندگار 

مدل HEC RAS‏ یا تحلیل سیستم رودخانه در مرکز مهندسی هیدرولوژی ساخته شد. این مدل توانایی 
تحلیل‌های هیدرولیکی جریان ۱ بعدی و ۲ بعدی در حالت‌های ماندگار و غیر ماندگار را دارد (رنگری" DATA‏ 
این پژوهش برای محاسبه پروفیل‌های سطح آب. پروفیل سرعت آب. طغیان سیلاب و عمق آب و محدوده پهنه‌های 
سیل گیر در دوره بازگشت‌های مختلف از تحلیل هیدرولیکی ۱ بعدی در حالت ماندگار در نرم‌افزار HEC RAS‏ 
نسخه T‏ استفاده شد. سپس با توجه به سیلاب‌های آنی در منطقه» شرایط مرزی برای سیستم رودخانه (بالادست و 
پایین‌دست) با توجه به عمق نرمال از نوع مختلط تعیین گردید. سپس مدل برای محاسبه پروفیل سطح آب در دوره 
بازگشت‌های de YO‏ ۱۰۰ و ۲۰۰ سال اجرا گردید. 

RAS mapper با ابزار‎ HEC RAS Jas تحلیل نتایج‎ —- Y-Y 

در نسخه‌های جدید مدل HEC RAS‏ ابزاری تحت عنوان RAS mapper‏ به مدل اضافه شده است. این ابزار 
توانایی تحلیل نتایج حاصل از مدل‌سازی ۱ بعدی و ۲ بعدی HEC RAS‏ برای هر دو حالت ماندگار و غیر ماندگار 


را دارد. ازاین‌رو در این پژوهش برای تحلیل نتایج (تعیین محدوده پهنه‌های سیل گیر و نقشه‌های az‏ سیلاب در 


1 Chow 
2 Rangari 


Lil a> YA‏ و مخاطرات محیطی شمارة سی و نهم 


دوره با ز گشت‌های ۵ Or‏ ۱۰۰ ۲۰۰ ساله) از ابزار RAS mapper‏ استفادہ گردید. در شکل 1 نمودار مراحل انجام 


این پژوهش آمده است. 


برآورد ded Bue uo‏ ایی و ا- برآورد دبی حداکٹر ded‏ با استفاده از ee‏ های متوسط حداکثر Shy‏ 

تحلیل فرلواقی میلاب 
Jd} Lei‏ ميلاب با Joli‏ توزیع های ehh‏ مناسب تردن توزیع با استقاده ‏ ازمون نکویی برازش در ترم IP‏ 
EasyFit‏ 


۱- ترسیم تشه توبوگرقی با امتفاده از نرم اقزار Global Mapper‏ و Google Earth‏ برای ترسیم e‏ رودخاته 


ترا ترس 


۲- ترسیم مجور اصلی رو دخاه و سولحل چپ و راست و عقاطع حرضی در (31 ۳11) AutoCAD‏ با استفاده از تقشه توبوگرافی 


۱- تقال عحور اصلی رودخاته سواحل جب و راست ترسیم شده در HEC RAS < AutoCADCivil 3D‏ 


۲- تعیین ضریب زبری Se‏ برای بتر رودخه و مواحل چپ و راست رو داه 
احرای HEC RAS Ju‏ ۳- وارد کرد دبی با دوره یازگشت های ۲۰۰۰۱۰۰۰۵۰۰۲۵ مه 


۴- تعین شرابط مرزی بالادست و بایین دست روتنخانه 


HEC RAS Ju با لجرای‎ jab تحلیل‎ -۵ 


۳- تعین پهته های سبل گر در دوره بازگشت های مختلف 


۴-ولید نقده خطر سیلاب در دورء بازگشت های۲۵ ۲۰۰۰۱۰۰۰۵۰۰ 
مه 

ux‏ حطر pls‏ روستای خنجه و باروتق S2‏ نوره 
com‏ ام abe‏ $22 


شکل -٤‏ روند نمای مراحل انجام این پژوهش 


۳- نتایج 
۱-۳- انتخاب بهترین توزیع و برآورد دبی حداکثر لحظه‌ای 
با توجه به نتایج حاصل جدول ۲ توزیع لوگ پیرسون نوع (LP3) Y‏ برای هر دو روش تجربی فیل - استینر و 
سنگال با توجه به آزمون کلموگروف اسمیرنوف برازش بهتری يافته است. در جدول Y‏ نتایج حاصل از آزمون نکوئی 
برازش جهت انتخاب بهترین توزیع آماری برای برآورد دبی حداکثر لحظه‌ای از دبی روزانه در نرم‌اف زار EasyFit5.6‏ 


OUS‏ داده شده است. هر چه مقدار آماره آزمون کلموگروف اسمیرنوف کمتر باشد نشان از برازش مناسب تابع 


سال دهم تهیه و تحلیل نقشه daz‏ سیلاب با استفاده از مدل .... Y4‏ 


توزیع نسبت به داده‌ها است (تان ' و همکاران:۲۰۲۱). نتایج حاصل از مقایسه مقدار ۲۷۵16 توزیع لوگ پیرسون 
نوع ۳ برای هر دو روش تجربی فیل - استینر و سنگال OUS‏ می‌دهد که روش تجربی روش فیل - استینر با بالاترین 
مقدار p-value‏ (معادل ۰/۹) مناسب‌ترین روش برای تبدیل دبی‌های حداکثر روزانه به دبی‌های حداکثر لحظه‌ای در 
ایستگاه هیدرومتری خنچه است. دلیران فیروز و همکاران (۱۳۹۶) نیز در مطالعه خود در حوضه آبخیز قمصر روش 
تجربی فیل - استینر را بهترین روش تجربی برای تبدیل دبی‌های حداکثر روزانه به دبی Sle‏ لحظه‌ای بیان کردند؛ 
بنابراین در این مطالعه برای تبدیل دبی‌های حداکثر روزانه به حداکثر لحظه‌ای از روش تجربی فیل - استینر و توزیع 
لوگ پیرسون نوع Y‏ استفاده شد. در جدول ۳ مقادیر برآورد شده دبی‌های حداکثر لحظه‌ای با روش تجربی فیل — 


استینر و توزیع لوگ پیرسون نوع ۳ آمده است. 


جدول ۲- نتایج حاصل برآورد دبی‌های حداکثر لحظه‌ای و تحلیل فراوانی سیلاب 


ردیف روش برآورد دبی لحظه‌ای توزیع آماری آماره آزمون کلمو گروف اسمیرنوف رتبه p- value‏ 

لوگ پیرسون نوع T Y‏ ۱ ۸ 

*/o£ Y JM پیرسون‎ 

rey Y Y گاما‎ 

"n í «£o نرمال‎ ze روش فیل‎ ) 
VIR! ۵ */00 Y پیرسون نوع‎ 

لوگ پیرسون نوع ۳ ۵ ۱ ۰/۸ 

eg Y 111 گاما‎ 

۲ Y -/í00 Y پیرسون نوع‎ 

"MT 1 */£04 روش سنگال نرمال‎ Y 
TEE ۵ ۰/۹ پیرسون‎ 


جدول ۳- نتایج برآورد دبی‌های حداکثر لحظه‌ای با روش فیل - استینر و توزیع لوگ پیرسون نوع Y‏ 


ردیف روش برآورد دبی حداکثر لحظه‌ای توزیع دوره باز گشت 
Ye UE os ۲۵‏ 
۱ فیل - استینر 1۳3 ۳/۱ WY ۱۹/۵ £v‏ 


] Tan 


" جغرافیا و مخاطرات محبطی شمارة سی و نهم 


HEC RAS شبیه‌سازی سیلاب در‎ yy 


در شکل ۵ پروفایل سطح آب در دوره باز £5 گشت Yo‏ ۵۰ ۱۰۰ ۲۰۰ ساله برای سه مقطع بالایی: میانی و پایینی 


رودخانه سوک چم به ترتیب در مقاطع عرضی شماره‌های ۰ و ۰ و ۰ شان داده شده | ست. بر اساس 


نتایج حاصل می‌توان نتیجه گرفت که از بالادست به سمت پایین‌دست شیب کم‌تر شده و بر سطح سیل گیری افزوده 


شده است. و پهنه‌های سیل uS‏ از دو طرف سواحل چپ و راست خارج می‌شود. 


T 
200 
Station (m) 


M 
WS TR-200 
WS TR=100 
سس‎ d - 
WS TR=50 
WS TR-25 
و‎ 
Ground 
Bank Sta 


Station (m) 


Bank Sta 


200 


Station (m) 


me .‏ آب برای سه مقطع (بالایی, میانی. پایینی) در دوره بازگشت (Vee ۷۰ ۰ YO)‏ ساله 


سال دهم تهیه و تحلیل نقشه daz‏ سیلاب با استفاده از مدل .... Y‏ 


oes -۳-۳‏ سطح سیل گیری در دوره با زگشت‌های مختلف 
با توجه به شکل 1 مساحت پهنه‌های سیل گیر در دوره بازگشست‌های مختلف ۲۵ (o‏ ۰ و ۰ سال به 
دوره باز گشت‌های مختلف YO)‏ 08( ۱۰۰ و (JL. Yes‏ در نرم‌افزار RAS mapper‏ آمده است. بر اساس eu‏ 


حاصل» از سراب به سمت پایاب رودخانه» با کاهش شیب و عریض شدن بستر بر مساحت سطح سیل گیری افزوده 


شده است. 


je‏ مساحت عرصه‌های سیل گیر با دوره با زگشت‌های YO‏ ۵۰, ۱۰۰ و ۲۰۰ ساله 


سطح سیل گبری برای دوره باژ گشت ۱۰۰ ساله 


سطح سبلل گیری برای عوره بازگشت ۲۰۰ اله 


شکل ۷- پهنه‌های سیل گیر در دوره بازگشت‌های Or YO)‏ ۱۰۰ و ۲۰۰ (all‏ در RAS mapper lo‏ 


۳۲ جغرافیا و مخاطرات محیطی شمارة سی و نهم 


۳-- تولید نقشه خطر سیلاب در دوره بازگشت‌های مختلف 

نقشه خطر سیل برای دوره بازگشت‌های ۲۵ ۵۰ ۱۰۰ ۲۰۰ ساله در محیط RAS mapper‏ برای بازه‌ای به طول 
۱ کیلومتر ترسیم شد؛ که نتایج آن در شکل‌های ۸ ٩‏ ۱۰ و ۱۱ نشان cols‏ شده است. مناطق در معرض خطر سیل 
در 7 كلاس Hi)‏ 48582 185 ) که به ترتیب نشان‌دهنده وضعیت امن امنیت بسیار کم کم» متوسط خطر زياد 
و خطر خیلی زیاد است تقسیم‌بندی شد (مک SU‏ ۰ ۲۰۱۶). در جدول ٤‏ طبقه‌بندی کلاس و میزان حطر سیلاب 
آمده است. همچنین در نشریه شماره ۳۰۷ راهنمای پهنه‌بندی سیل و تعیین حد بستر و حریم رودخانه دوره بازگشت 
ale ۰‏ به‌عنوان پهنه سیل پرخطر معرفی گردید. 

نتایج حاصل از تهیه نقشه‌های خطر به دست آمده در این پژوهش OUS‏ می‌دهد که در تمام دوره بازگشت‌ها 
«en‏ دو روستای خنچه و بارونق در معرض سیل قرارگرفته و در کلاس‌بندی متوسط تا خیلی زیاد (H6-HA)‏ (ناامن 
برای وسیله نقلیه و مردم و ساختمان‌ها و زمین‌های در معرض خطر سیلاب) قرار دارند. لذا اقدامات آبخیزداری و 
کنترل سیل در این دو روستا ضروری است. همچنین بین سال‌های ۱۳۸۰-۱۳۳۰ در شهرستان کاشان تعداد V4‏ 
سیلاب رخ داده Cul‏ که خسارت‌های زیادی به همراه داشته است (معیری و انتظاری» ۱۳۸۷). ازآنجایی که هیچ گونه 
عملیات کنترل سیلاب و آبخیزداری در جهت مهار سيل برای این دو روستا صورت نگرفته» لذا از کل مساحت دو 
روستای خنچه و بارونق به ترتیب ۹/۸۸ و ۱۶/۳ هکتار حدود 1/۱ و ۰/۱ هکتار به ترتیب 1۱/۵ و 1/۵ درصد از 


اراضی این دو روستا در معرض خطر سیلاب با دوره yb‏ گشت‌های ۰ ساله )» Ce‏ قرار دارند. 


حدول -٤‏ طبقه‌بندی کلاس و میزان خطر سیلاب (مک لاک ی" QM‏ 


ردیف کلاس خطر میزان خطر توصیف 
Hı ۱‏ امن به‌طور کل ایمن برای مردم و وسیله نقلیه و ساختمان‌ها 
H2 Y‏ بسیار کم ناامن برای وسیله نقلیه کوچک 
e$ H3 ۷‏ ناامن برای وسیله نقلیه و بچه‌ها و افراد سالخورده 
H4 t‏ متوسط ناامن برای وسیله نقلیه و مردم 
E Hs ۵‏ ناامن برای وسیله نقلیه و مردم و بعضی از ساختمان‌های در معرض خطر 
He E‏ خیلی زیاد ناامن برای وسیله نقلیه و مردم و همه ساختمان‌های در معرض خطر 


] McLuckie 
2 McLuckie 


۳ .... سیلاب با استفاده از مدل‎ daz و تحلیل نقشه‎ ang 
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شکل 4- نقشه خطر سیل برای دوره بازگشت *0 ساله در محیط RAS mapper‏ 


۳۶ 
جغرافیا و مخاطرات محیطی 
شمارة سی و نهم 


۱- نقشه > s‏ 
شکل نقشه خطر سیل برای دوره باز گشت ۲۰۰ ساله در محیط RAS mapper‏ 


سال دهم تهیه و تحلیل نقشه daz‏ سیلاب با استفاده از مدل .... Yo‏ 


ci‏ جمع‌بندی 

در این پژوهش برای تولید نقشه‌های خطر سیل در دوره بازگشت‌های مختلف ابتدا نقشه توپوگرافی منطقه با 
استفاده از نرم‌افزار google earth‏ و Global mapper‏ تهیه گردید. سپس برای ترسیم مقاطع عرضی» پروفیل طولی 
رودخانه و ایجاد سطح از منطقه نقشه توپوگرافی تهیه شده به نرم‌افزار AutoCAD Civil 3D‏ انتقال داده شد. سپس 
مقاطع عرضی ساخته شده به نرم‌افزار HEC RAS‏ انتقال داده شد. همچنین ازآنجایی که منطقه موردمطالعه در منطقه 
خشک واقع شده است و فاقد آمار دبی‌های حداکثر لحظه‌ای است و تنها دارای آمار دبی‌های حداکثر روزانه است 
برای تبدیل دبی‌های حداکثر روزانه به حداکثر لحظه‌ای با توجه به مقدار آماره آزمون کلموگروف اسمیرونوف 
(۰/۱۳) و مقدار p-value‏ (۰/۹) می‌توان نتیجه گرفت روش فیل - استینر روشن مناسبی است. سپس آنالیز فراوانی 
سیلاب با استفاده از ۵ توزیع (لوگ پیرسون نوع Y‏ پیرسون. US‏ (حداکثر» نرمال» پیرسون نوع ۳) در نرم‌افزار 
EasyFit‏ مورد تحلیل قرار گرفت. 

با توجه به نتایج آنالیز فراوانی سیلاب از توزیع لوگ پیرسون نوع ۳ برای ایجاد دبی در دوره بازگشت‌های Y0)‏ 
(Yer cher ۰‏ ساله استفاده گردید. سپس با استفاده از نتایج شبیه‌سازی سیلاب در دوره بازگشت‌های فوق در RAS‏ 
۲ نقشه‌های خطر سیل در دوره بازگشت‌های ۲۵ ۰٠۰۰ Or‏ ۲۰۰ ساله برای بازه مطالعاتی به طول ۱۱ کیلومتر 
رودخانه سوک چم تهیه گردید. از آنجایی در نشریه شماره ۳۰۷ برای پهنه‌بندی سیل و تعیین حد بستر رودخانه 
دوره باز گشت ۱۰۰ ساله را به‌عنوان پهنه‌های سیل پرخطر معرفی شده است لذا در این پژوهش برای تحلیل عناصر 
در معرض خطر سیلاب دو روستای (خنچه و بارونق) برای دوره بازگشت ۱۰۰ ساله محاسبه گردید و برای سایر 
دوره بازگشت‌ها (۲۵ ۵۰ ۲۰۰) ساله تنها نقشه‌های حطر تولید گردید. ازآنجایی‌که دو روستای خنچه و بارونق در 
امتداد بازه مطالعاتی قرار دارند و هیچ‌گونه عملیات کنترل سیلاب و آبخیزداری برای این دو روستا صورت نگرفته و 
al pos‏ در تمامی دوره بازگشت‌ها در معرض خطر سیل قرار دارند بااین‌وجود از نتایج تحلیل خطر سیلاب برای دوره 
بازگشت ۰ ساله Gha p)‏ برای این دو روستا نشان از این دارد که از کل مساحت دو روستای خنچه و بارونق به 
ترتیب با ۹/۸۸ و ۱۶/۱۳ هکتار حدود 2/۱ و ۰/۱ هکتار در معرض سیلاب قرار می گیرند و در کلاس خطر سیلاب 
متوسط تا خیلی زیاد (ناامن برای وسیله نقلیه و مردم و ساختمان‌های در معرض خطر) قرار دارند. همچنین مساحت 
سیل گیری برای دوره باز گشت‌های ۲۵ ۵۰ ۰۱۰۰ ۲۰۰ ساله به ترتیب ۰/۷۷ ۰/۹ ۰۱/۱۵ ۱/4۷ کیلومترمربع به دست 
آمده است. لذا ایجاد دیواره‌های ساحلی و سامان‌دهی رودخانه در طول مسیر عبور رودخانه از این دو روستا امری 
ضروری است که می‌تواند خطرات ناشی از سیل را تعدیل نماید. استفاده از نقشه‌های تهیه شده برای انجام اقدامات 
کنترل سیل» سامانه‌های پیش‌بینی و هشدار سیل بیمه سیل و به‌طورکلی جهت اقدامات سازه‌ای و غیرسازه‌ای 
مهار سیل در منطقه پيشنهاد می‌شود. 
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